Dissipation induced Tonks-Girardeau gas of photons 
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A scheme for the generation of a Tonks-Girardeau (TG) gas of photons with purely dissipative in- 
teraction is described. We put forward a master equation approach for the description of stationary 
light in atomic four-level media and show that, under suitable conditions, two particle decays are 
the dominant photon loss mechanism. These dissipative two-photon losses increase the interaction 
strength by at least one order of magnitude as compared to dispersive two-photon processes and 
can drive the photons into the TG regime. Our scheme allows for measurements of various char- 
acteristic correlations of the TG gas via standard quantum optical techniques, including quantities 
that distinguish it from free fermions. 
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Quantum mechanics categorizes particles into fermions 
or bosons. In three dimensions only these two categories 
are possible, whereas more exotic anyons can exist in two 
dimensions In one dimension, the particle statistics 
can not be considered without taking inter-particle in- 
teractions into account 0- A prominent example arc 
bosons that interact via strong repulsive forces in a one- 
dimensional setting and can enter a Tonks-Girardeau 
(TG) gas regime [3|, where they behave with respect to 
many observables as if they were fermions. A TG gas 
can be described as the strong interaction limit of the 
Lieb-Liniger model Q]. 

Strong correlations in many-particle systems, such as 
in the TG gas, give rise to interesting and partly not 
yet well understood physics. A substantial amount 
of research is thus currently devoted to these systems 
and progress in cooling and trapping of atoms and ions 
has opened up possibilities to study strongly interacting 
many-body systems experimentally with unprecedented 
precision. Eventually, this progress enabled the obser- 
vation of a TG gas of atoms in an optical lattice Q. 
Later, an experiment with cold molecules showed that 
not only elastic interactions, but even two-particle losses 
alone are able to create a TG gas where two molecules 
never occupy the same position, thereby avoiding dissipa- 
tion of particles. This counterintuitive result can be re- 
garded as a manifestation of the quantum Zeno effect Q . 

In contrast to atoms, photons are massless particles 
that do not interact at all. Nonetheless, effective many- 
body systems of photons and polaritons can be generated 
by employing light matter interactions. This concept has 
been introduced recently [ol . [lol. 
receiving increasing attention |12l . 
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and is currently 
To enter 



the strongly correlated regime and access its rich physics, 
sufficiently strong effective interactions are needed. Suit- 
able experimental setups, like arrays of coupled microcav- 
ities doped with emitters [l3| or optical fibers that couple 
to atoms 1^ |l7[ , thus need to combine strong photon- 
emitter coupling and low-loss photon propagation. In all 
these setups, the main challenge for realizing strong cor- 
relations is to make the polariton-polariton interactions 



much stronger than photon losses which are inevitably 
present in every experiment. 

Here we present an effective many-body system of 
polaritons where the ubiquitous but usually undesircd 
dissipative processes become the essential ingredient for 
the creation of strong many-particle correlations. This 
paradigm shift allows us to relax some conditions on the 
model parameters such that the achievable nonlincarities 
in our approach arc at least an order of rnagnitude larger 
than their conservative counterparts [1, ij, 14 1. In par- 
ticular, we show that the dissipative nonlincarities in our 
system give rise to a TG gas of photons. For this regime, 
fermionic (e.g. Friedcl oscillations) as well as non-local 
(e.g. the single particle density matrix) correlations of 
the TG gas can be measured via standard quantum op- 
tical techniques. 

We consider photons guided in an optical fiber that 
interact with nearby atoms [l^ . 17 1, where stationary 
light [1^ is created via Electromagnetically Induced 
Transparency (EIT) [l^. As compared to coupled mi- 
crocavities, the fiber approach is appealing due to the 
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FIG. 1: (Color online) (a) Considered setup of N atoms con- 
fined to an interaction volume of length L and transverse 
area A. Q,± axe the Rabi frequencies of the classical control 
fields, and E± axe the quantum probe fields, (b) Atomic level 
scheme. 7ij is the full decay rate on the \i) <-> \j) transition, 
5 and A label the detuning of the probe fields with states |3) 
and |4), respectively, and e is the two-photon detuning. 
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low photon loss of the fiber and since the longitudinal 
trapping of light is done optically, thus avoiding the need 
to build many mutually resonant cavities. We thus focus 
on this setup here. However, our mechanism for building 
up correlations works equally well in cavity arrays, and 
the dissipative nonlinearities we discuss here are always 
stronger than their conservative counterparts indepen- 
dent of the geometry of the experimental device [20|. 

We start with a more detailed description of our one- 
dimensional model shown in Fig.[T] Each of the N atoms 
interacts with control and probe fields denoted by f2± and 
E±, respectively. The control fields of frequency ujc are 
treated classically and fl^ {^-) labels the Rabi frequency 
of the control field propagating in the positive (negative) 
z direction. In addition, we assume that the control fields 
are spatially homogeneous but may depend on time. The 
probe fields E+ and E- are quantum fields that propa- 
gate in the positive and negative z direction, respectively. 
They are defined as E±{z) = J2k '^iKe'^^^^ , where a±K 
are photon annihilation operators. The wave numbers K 
are positive and of the order of the wave number kc of 
the control field. 

We model the time evolution of the atoms and the 



quantized probe fields by a master equation [2lJ for their 
density operator g, g = —jr[H, g] + C^g, where C^g de- 
scribes spontaneous emission from states |3) and |4), and 
the full decay rate on the transition \i) <-> |j) is denoted 
by 7y (see Fig. [T|) . In a rotating frame that removes the 
time-dependence of the classical laser fields, the system 
Hamiltonian H reads H = Hq + H\ + -f^NL, where 

Ho = - h^{ujp - luk) (a^ax + a^A-o-A' 
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(1) 



describes the free time evolution of the atoms and the 



probe fields. A^^ 



i){ifi\ is a projection operator onto 



state of atom n, the energy of level |i) is huji (we 
set LOi =0), and transition frequencies are denoted by 
We denote the central frequency of the 
probe pulse by LUp. The detuning of the probe field with 
respect to the transitions |3) ^ |1) and |4) ^ |2) is 



1 and A 



W42, respectively, 



labeled by S — l 
and e — {tOp ~ ojc) ~ ^2 is the two-photon detuning. The 
interaction between the atoms and the probe and control 
fields is described by Ha + -ffNL, with 
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NL 



+ giSi^^ [s+iz,) + E^{z^)] } + h.c. , (2) 
E+{Zf,)+E^{Zf,)\+h.c.. (3) 
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Transition operators of atom fi at position are de- 



fined as S, 



Iv)Oa'I (* 7^ i)' and gi and 52 are 



the single-photon Rabi frequencies on the |3) <-> |1) and 
1 4) ^ 1 2) transitions, respectively. In the following, we 
assume that the Rabi frequencies of the control fields 
are identical (and real) and set r2+ = r2_ = ilc- With 
this choice, the interaction of the probe and control fields 
with the A-subsystem formed by states |2) and |3) al- 
lows to store the probe field inside the medium [iSj. On 
the other hand, the coupling of the probe fields to the 
1 4) ^ 1 2) transition creates an effective photon-photon 
interaction [il ]. 

Next we outline the approach we developed to reduce 
the master equation g = —j-[H, g] + C^g for the atoms 
and quantized probe fields into a master equation solely 
for dark-state polaritons [2^, formed by collective exci- 
tations of photons and atoms. We represent the quan- 
tum state of dark-state polaritons by a density matrix gD 
comprised of dark-states \a) = Y\k=ii^ / V^iJ-) {i/jI)"''' \0) 
that satisfy HaIch) = 0, and the vacuum state |0) = 
|{0}phot; li, . . . , Iw) is the state where all photon modes 
of the probe fields are empty and all atoms are in state 
The operators ipk are defined as 24 1 
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(4) 



where sin 6 = 



Wgi/ilo, COS 6* = ^/2Vtc/^o^ and ilo = 
The operator Ak = (a^^+fe -I- a-k,+k)/V^ 



is a superposition of two counterpropagating probe field 
modes, and describes the spin coherence, = 
-ikzu . Note that the wave number k can 
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be positive or negative, and for all relevant k we have 
|fc| ^ kc- We assume that initially all atoms are in state 
|1) and that the total number of photons is much smaller 
than the number of atoms N. In this case, the dynamics 
induced by the Hamiltonian H is confined to a subspace 
Hfe of the total state space where {-tpl ^11^ IV") ~ ^ 

for all G TipE- It follows that the operators tpk obey 
bosonic commutation relations in TipE, = ^kp, 

where we neglected corrections of order 1/N. 

The dark-state polaritons are eigenstates of Ha, but 
the remaining parts Hq and i?NL of the system Hamilto- 
nian give rise to a non-trivial time evolution of gu . Fortu- 
nately, this dynamics can be studied entirely in terms of 
bosonic quasi-particle excitations if the system dynamics 
is restricted to the subspace TYfe- In particular, the free 
time evolution Hq introduces a coupling of dark-state po- 
laritons to bright polaritons, (j)k = ^fcsin^ -I- Xfacos^, 
and photons, = {ak^+k ~ O'-k^+k) I ^f^- These excita- 
tions are in turn coupled to the excited state |3). Fur- 
thermore, i?NL introduces a direct coupling of dark-state 
polaritons V'fc to the excited state |4) via a two-particle 
process [2^. Excitations in the states |3) and |4) are ere- 



ated by p2,+ , 



and f/fct ^. 



± 



U\ _ , respectively, where 

i(ka-k)z^^j^2'M 



3 



t 



Finally, we note that spontaneous emission from states 
1 3) and |4) results in the decay of excitations created by 
Pl± and Ul^. 



We employ projection operator techniques [2l| to de- 
rive a master equation for gD which is obtained from g 
by a partial trace over all excitations except for the dark 
state polaritons -^fe. We restrict our analysis to the so- 
called slow-light regime where sin^ « 1 and cos^ <SC 1. 
In this case, the coupling of dark-state polaritons to ex- 
citations in state |3) and |4) is much slower than the 
decay of the relevant correlation functions (0fc0p(T), 
(Z)fcZ)t)(r), and (Uk,+Up^^){T), which happens on a 
timescale given by the lifetimes of the exited states |3) 
and |4). This existence of two different time scales allows 
us to derive a master equation in Born-Markov approx- 
imation if 4.91 cos^ 6'iVph < 7I2, cos'^ 9 c^k"^^ J nl < 1, 
cos^ OAuj'^ /^Iq <^ 1, and ^Iq ^ 7^ , \5\. Here c is the speed 
of light, Alu ~ LOp^LOc is the frequency difference between 
the probe and control fields and A'ph is the number of 
photons in the pulse. We describe the polariton pulse by 
the field operator = X]fe e^'^^V'fc which obeys 

the commutation relations [-(/'(z), '(/'^(^Ol = 5{z — z'). The 
maximal wave number contributing to ip is fcmax- Fur- 
thermore, our derivation assumes that fast oscillating 
spin coherences with wave number ±2A:c are washed out 
due to atomic motion [18j . For a small two-photon de- 



tuning £ 



cos 9Alo, wc obtain 



hgx^ = -iHcsgo + igiyHlff + Ig^, + Cig-D+ Hig-D-, (5) 
where _ffcff is a non-hermitian Haniiltonian, 



off : 



2m, 



off Jo 



^dza,Vt5,^ + |^''dz^t2^2^ (6) 



TOoff = —hVi^ I (^bt? co^ Q) is the effective mass of the 
polaritons, g — 2h,Lg\ cos^ ^ ~ cos^ dAto + 1742/2) is 
the complex coupling constant, and 



Tgxy = -Im(.g) / dzV £'dV' 
Jo 

hTAuj^Vii;] ^ 



(8) 



Here V[X] = /^^ Az{X^XgY) + gnX'^ X - 2XgjiX^) is a 
dissipator in Lindblad form [2l| for an operator X, and 
r = 731 -I- 732 is the full decay rate of state |3). For 
optical fibers, photon losses due to leakage are very low 
and can be neglected. If they need to be taken into ac- 
count, an additional decay term with the same structure 
as vC2^?D but with a decay rate kcos^ 9 appears (k is the 
bare photon leakage rate). To confirm the accuracy of 
our results, we compared the predictions of the master 
equation (O for the A- subsystem {g = 0) to the results of 
a full numerical integration of MaxwcU-Bloch equations 
for classical fields and found excellent agreement. 



Next we derive the essential results of this letter from 
the master equation ([5]) that describes a one-dimensional 
system of interacting bosons. The first contribution to 

TJeff in Eq. ([6]), {h^/2mes) dzdzij^^dzip, represents a 
kinetic energy term with quadratic dispersion relation 
for the polaritons. The term proportional to g in Eq. ([5]) 
and in Eq. Q account for elastic and inelastic two- 
particle interactions that originate from the coupling of 
dark-state polaritons to the excited state |4). More pre- 
cisely, the real part of g gives rise to a hermitian con- 
tribution to HeS that accounts for elastic two-particle 
collisions. On the other hand, the imaginary part of g 
together with XpD gives rise to a two-particle loss term 
that can be written in Lindblad form as —lm{g/2)'D[ip'^]. 

The contributions Cigu and £2^0 describe single- 
polariton losses that can be omitted under the following 
conditions. Since Cigi^ is proportional to Aw^, single- 
particle losses are minimized by minimizing |Aa;|. Note 
that this fact has not been pointed out so far. From now 
on we assume that Auj'^ is small enough such that C2gu 
represents the dominant single particle losses. This is 
reasonable if \Auj\ is at most of the order of GHz and 
implies \e\ <C |724|- The term £2^0 is negligible if two 
conditions arc met. First, the dynamics induced by the 
kinetic energy term proportional to meff in Eq. ([6]) must 
be fast as compared to the inverse decay rate of polari- 
tons introduced by £290- This can be achieved if we set 
\5\ 3> r. Second, losses due to £2^0 must be negligi- 
ble which imposes a limit on the maximal evolution time 
^max *C 2fll/{Tc'^k'^^^ cos^ 9). This implies that fmax can 
be of the order of l/(cos2 9r) > l/F. 

Under these conditions, the master equation ([S]) re- 
duces to ftpb = —iHcsQ'D + iQuHlff + ^Qu and can be 
identified with the generalized Lieb-Liniger model 0] for 
a one-dimensional system of bosons with mass Woff and 
complex interaction parameter g. All features of the 
Lieb-Liniger model 0, 0] are characterized by a single, 
dimensionless parameter G = 'mcsg/{h^Np\i/L), where 
iVph is the number of photons in the pulse. The absolute 
value of G is 



\G\ 



(l/16)r7420D2 
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where OD = ANglL/{cT) = 4iV5|L/(c742) is the op- 
tical depth on the probe field transitions. Note that 
the parameters giL and gfL are independent of the 
length of the system since gi , (72 ~ 1 / VAL. It fol- 
lows that the parameter G and the optical depth de- 
pend only on the transverse area A of the interaction 
volume, but not on the length L of the cell. The 
absolute value of G characterizes the effective interac- 
tion strength between the particles. In the strongly 
correlated regime \G\ ^ 1, the interaction between 
the particles creates a Tonks-Girardeau gas where pho- 
tons behave like impenetrable hard-core particles that 
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never occupy the same position. Formally, this re- 
sult can be derived via the pair correlation function 
5(2)(z,z') = {i^^{z)i;^{z')i,{z)i,{z'))/{{h{z)){n{z'))) with 
h{z) = il}\z)-4>{z). For the ground state of the general- 
ized Lieb-Liniger model in the strongly correlated regime, 
5(2) (2, z) = (1 - l/^ph)47rV(3lG|2) is close to zero and 
vanishes in the limit |G| — > 00 [7[. Moreover, this ground 
state is the same as in the original model with repul- 
sive interaction for \G\ 00. It follows that g'^'^\z,z') 
for z ^ z' exhibits Friedel oscillations [2^ that indicate 
a crystallization of photons in the fiber. 

The parameter \G\ is maximal if the interaction be- 
tween the polaritons is purely dissipative (A = 0). Since 
the realization of a regime where the two-particle in- 
teractions are dominated by elastic processes requires 
A ^ 742/2, the conservative nonlinearities are at least 
an order of magnitude smaller than the dissipative coun- 
terparts for A = 0. It follows that purely dissipative 
interactions between the polaritons we discuss here are 
most effective for the generation of correlations. 

An analysis of dissipation-induced correlations (A = 
0) requires at least two photons. Assuming \5\/T = 10 
such that the single-particle loss term £igD in Eq. ([5]) is 
negligible and iVph = 2, Eq. © shows that \G\ is larger 



17| 



than unity for OD^/iV > 160. A recent experiment 
with atoms loaded into a hollow fiber reports a value of 
OD^/N « 0.3. If we assume for simplicity that the decay 
rates of the atomic states |3) and |4) do not depend on 
the transverse area A, we find IG] oc NXp/A"^, where Ap is 
the wavelength of the probe field. It follows that \G\ does 
not dependent on the strength of the atomic transition 
dipole moments and could be increased by a reduction of 
the area A or by an increased number of atoms N inside 
the fiber. In contrast to cavity QED systems [Ol, we 
point out that the condition (72 3> 742 is not required to 
obtain large values of \G\. 

The observation of the dissipation-induced TG gas 
regime requires that the system can be prepared in 
low-energy states. One possibility is the procedure 
described in [l2l which relies on an adiabatic state 
transfer realized by a time-dependent detuning A. A 
second possibility does not require any tuning of the 
two-particle losses. The master equation ([Sj implies 
that losses due to inelastic two-particle interactions 
are related to the pair-correlation function g^^'>{z,z) 
via dt{n{z)) = {2/h)lm{g)g^^\z,z){n{z)f . It follows 
that uncorrelated states with g^'^^z,z) w 1 decay 
much faster than those where g'-'^\z,z) w 0. There- 
fore, a regime where g^^\z,z) < 1 should be en- 
tered on a time scale h/[21m{g)Np\i/ L], which is shorter 
than the maximally allowed evolution time tmax if 
165r|7Vphf2o/(r742C^^max) > 1- Siucc the ground state of 
the generalized Lieb-Liniger model decays at the smallest 
rate 0], two-particle losses themselves are then able to 
drive the system into states close to the ground state. A 
more rigorous investigation of this point would require a 



numerical integration of the master equation ([5]) which 
is beyond the scope of this work. 

For measurements, we note that the polariton pulse 
can be released from the fiber without distortion if 
one control field is adiabatically switched off 
It follows that spatial correlations {ip'^ {z)ip{z')) and 
(ip^ {z)ip^z')ip{z)^p{z')) of the trapped pulse are mapped 
into first and second order correlations in time of the out- 
put light, respectively. Since the latter can be detected 
via standard quantum optical techniques, Friedel oscilla- 
tions of g^^^Zjz'), the correlations {'ilj^z)ip{z')) and the 
characteristic momentum distribution of the TG gas, can 
be measured with high precision. 

The authors thank W. Zwerger for discussions. This 
work is part of the Emmy Noether project HA 5593/1-1 
funded by the German Research Foundation (DFG). 
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